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The  effect  of  the  cathode  catalyst  layer’s  structure  and  composition  on  the  overall  performance  of  a  poly¬ 
mer  electrolyte  membrane  fuel  cell  (PEMFC)  is  investigated  numerically.  The  starting  point  of  the  sub-grid 
scale  catalyst  layer  model  is  the  well-known  flooded  agglomerate  concept.  The  proposed  model  addresses 
the  effects  of  ionomer  (Nation)  loading,  catalyst  (platinum)  loading,  platinum/carbon  ratio,  agglomerate 
size  and  cathode  layer  thickness.  The  sub-grid  scale  model  is  first  validated  against  experimental  data  and 
previously  published  results,  and  then  embedded  within  a  two-dimensional  validated  computational 
fluid  dynamics  code  that  can  predict  the  overall  performance  of  the  fuel  cell.  The  integrated  model  is  then 
used  to  explore  a  wide  range  of  the  compositional  and  structural  parameter  space,  mentioned  earlier.  In 
each  case,  the  model  is  able  to  correctly  predict  the  trends  observed  by  past  experimental  studies.  It  is 
found  that  the  performance  trends  are  often  different  at  intermediate  versus  high  current  densities— the 
former  being  governed  by  agglomerate-scale  (or  local)  losses,  while  the  latter  is  governed  by  catalyst  layer 
thickness-scale  (or  global)  losses.  The  presence  of  an  optimal  performance  with  varying  Nation  content 
in  the  cathode  is  more  due  to  the  local  agglomerate-scale  mass  transport  and  conductivity  losses  in  the 
polymer  coating  around  the  agglomerates  than  due  to  the  amount  of  Nation  within  the  agglomerate.  It 
is  also  found  that  platinum  mass  loading  needs  to  be  at  a  moderate  level  in  order  to  optimize  fuel  cell 
performance,  even  if  cost  is  to  be  disregarded. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

One  of  the  impediments  to  commercial  success  of  fuel  cells  is 
high  cost.  One  of  the  major  contributing  factors  to  high  cost  is  the 
excessive  use  of  platinum  within  the  catalyst  layers.  Thus,  there  is 
strong  motivation  for  using  as  little  platinum  as  possible  without 
hampering  the  overall  performance  of  the  fuel  cell.  In  a  polymer 
electrolyte  fuel  cell  (PEMFC),  since  the  catalyst  (typically  platinum) 
is  dispersed  within  a  complex  porous  matrix  comprised  of  carbon, 
the  ionomer  (typically  Nation)  and  platinum,  the  overall  perfor¬ 
mance  of  the  fuel  cell  does  not  scale  linearly  with  an  increase  in 
the  amount  of  platinum  within  the  catalyst  layer.  Rather,  the  per¬ 
formance  depends  on  how  effectively  the  platinum  is  dispersed 
within  the  porous  matrix  so  that  most  of  it  is  actually  utilized  in  cat¬ 
alyzing  the  electrochemical  reactions.  Since  the  proper  functioning 
of  a  PEMFC  cathode  requires  existence  of  triple  phase  bound¬ 
aries  [1,2]  between  the  Nation  (for  proton  transfer),  platinum  (for 
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catalysis)  and  carbon  (for  electron  transfer),  the  determination  of 
optimum  composition  and  structure  of  the  catalyst  layer  is  a  monu¬ 
mental  task.  Such  studies  are  often  undertaken  using  experiments 
[3,4],  but  are  very  time-consuming  and  expensive.  Needless  to  say, 
with  trial-and-error  experimental  procedures,  only  a  small  portion 
of  the  overall  design  space  can  be  explored  within  a  reasonable 
time.  Computational  modeling  provides  an  alternative  means  to 
address  this  critical  issue.  However,  such  models  also  require  care¬ 
ful  formulation  and,  to  some  extent,  calibration  of  unknown  fitting 
parameters.  Nevertheless,  given  the  importance  of  the  problem  in 
the  context  of  fuel  cell  technology  advancement,  and  the  difficul¬ 
ties  associated  with  experimental  trial-and-error  procedures,  there 
is  strong  motivation  to  try  to  address  the  problem  using  models, 
either  fundamental  or  phenomenological. 

Over  the  past  decade  or  so,  there  has  been  an  explosion  in  the 
number  of  technical  publications  that  have  reported  findings  of  the 
simulation  of  the  catalyst  layer  of  a  PEMFC.  Broadly,  these  simula¬ 
tions  may  be  classified  into  two  categories:  (1)  simulations  that 
utilize  approximate  models  based  on  the  pseudo-homogeneous 
film  concept  [5-9],  or  the  flooded  agglomerate  concept  [10-24], 
and  (2)  direct  numerical  simulation  (DNS)  of  the  catalyst  layer 
[25,26].  While  the  latter  approach  is  more  powerful,  it  has  seri¬ 
ous  limitations.  First,  it  requires  reconstruction  of  the  catalyst  layer 
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Nomenclature 

Av  total  catalyst  surface  area  per  unit  volume  of  cath¬ 
ode  (m-1) 

A0  total  catalyst  surface  area  per  unit  mass  of  catalyst 
particle  (m2  kg-1 ) 

c  dissolved  oxygen  concentration  in  Nation 

(kmolm-3) 

c*  dissolved  oxygen  concentration  in  Nation  in  equilib¬ 

rium  with  inlet  gas  (kmol  m-3) 

Co2,g  oxygen  gas  concentration  in  cathode  gas  pores 
(kmolm-3) 

cq2  g  oxygen  gas  concentration  at  cathode  inlet 
(kmolm-3) 

cjjef  standard  reference  oxygen  concentration 

(kmolm-3) 

Cw  concentration  of  water  in  membrane  (kmol  m-3 ) 

d  average  pore  size  of  cathode  (m) 

Dkn  binary  diffusion  coefficient  of  species  k  into  n 

(m2  s-1) 

Do2,n  diffusion  coefficient  of  oxygen  in  Nation  (m2  s-1 ) 
DT  temperature  dependence  of  (m2  s-1 ) 

Dw  diffusion  coefficient  of  water  (m2  s-1 ) 

Dx  diffusion  coefficient  of  water  (m2  s-1 ) 

D'x  concentration  dependence  of 

Dol  n  effective  diffusion  coefficient  of  oxygen  in  Nation  in 
agglomerate  (m2  s-1 ) 

Effi  nucleus  effectiveness  factor 

Eff2  film  effectiveness  factor 

F  Faraday’s  constant  (96.487  x  106  C  kmol-1 ) 

i  net  current  density  vector  (A  m-2 ) 

ip  ionic  phase  current  density  vector  (A  m-2 ) 

is  electronic  phase  current  density  vector  (A  m-2 ) 

z£at  surface  current  density  on  a  smooth  catalyst  surface 

(Am-2) 

igef  standard  exchange  current  density  on  cathode 

(Am-3) 

jo  reference  current  density  (A  m-2 ) 

jf1  net  transfer  current  at  anode  due  to  electrochemical 

reaction  (Am-2) 

j£at  net  transfer  current  density  at  cathode  due  to  elec¬ 
trochemical  reaction  (Am-3) 

Jadv  water  flux  due  to  advective  transport 

(kmolm-2  s-1) 

Jdiff  water  flux  due  to  diffusion  (kmol  m-2  s-1 ) 

Jdrag  water  flux  due  to  electro-osmotic  drag 

(kmolm-2  s-1) 

]k  diffusion  mass  flux  of  the  kth  species  (kg  m-2  s-1 ) 

L  cathode  catalyst  layer  thickness  (m) 

mpt  platinum  mass  loading  (kg  m-2 ) 

Mm  molecular  weight  of  the  membrane  (kg  kmol-1 ) 

Mk  molecular  weight  of  kth  species  (kg  kmol-1 ) 
n  number  of  electrons  transferred  during  the  electro¬ 

chemical  reaction 

h  number  of  agglomerates  per  unit  volume  of  cathode 

(m-3) 

N  total  number  of  gas-phase  species 

p  pressure  (Pa) 

Pt|C  platinum-carbon  mass  ratio  in  catalyst  layer  ink 
ragg  radius  of  agglomerate  (m) 

R  universal  gas  constant  (8314J  kmol-1  K-1 ) 

Sk  production  rate  of  kth  species  (kg  m-3  s-1 ) 


(S/V)e ff  Effective  surface  to  volume  ratio  (m-1 ) 

t  time  (s) 

T  temperature  (K) 

U  Bulk  fluid  velocity  (ms-1) 

V0c  open  circuit  voltage  (V) 

Yk  mass  fraction  of  kth  species 

Greek  letters 

ota,ac  Tafel  constants  for  anode 

aT  Tafel  constant  for  cathode  catalyst  model 
Pk  concentration  exponents  for  the  kth  species 
8  polymer  coating  thickness  around  agglomerate 

nucleus  (m) 
s  wet  porosity 

£agg  volume  fraction  of  polymer  in  agglomerate  nucleus 
£cat  porosity  of  cathode  catalyst  layer 

£S  volume  fraction  of  platinum  +  carbon  in  cathode 

£N  volume  fraction  of  polymer  in  cathode 

rj  electrode  overpotential  (V) 

p d  electro-osmotic  drag  coefficient 

k  permeability  ( m2 ) 

X  water  content 

[Ak\  molar  concentration  of  species  k  (kmolm-3) 

li  dynamic  viscosity  (kg  m-1  s-1 ) 

p  mass  density  of  mixture  (kg  m-3 ) 

pc  density  of  carbon  (kg  m-3 ) 

pPt  density  of  platinum  (kg  m-3 ) 

density  of  dry  membrane  (kg  m-3 ) 
a  electrical  conductivity  (£2-1  m-1 ) 

(7 p  electrical  conductivity  of  the  ionic  phase  ( £2-1  m-1 ) 

cr c  electrical  conductivity  of  the  electronic  phase 
(^-1m-1) 

<t30  concentration  dependence  of  electrical  conductivity 

(ft-1  m-1) 

rcat  tortuosity  of  cathode 

<2>f  ionic  phase  potential  (V) 

0S  electronic  phase  potential  (V) 


microstructure  from  SEM  or  TEM  images.  Secondly,  direct  numer¬ 
ical  simulations  are  difficult  to  perform  and  are  very  expensive. 
Simulation  times  for  a  single  case  may  often  run  into  days.  Finally, 
certain  compromises  have  to  be  made  with  regard  to  the  actual 
geometry  [25],  rendering  the  results  of  these  studies  inconclu¬ 
sive  in  many  cases.  Nevertheless,  this  approach  is  often  useful 
for  fundamental  understanding  of  the  coupling  between  transport 
phenomena  and  heterogeneous  reactions  at  the  micro-scale  and 
can  be  used  as  a  framework  to  develop  simplified  models. 

The  more  popular  approach  is  based  on  models  of  coupled  trans¬ 
port  and  reactions  within  the  catalyst  layer  structure.  Historically, 
two  different  model  types  have  been  used  for  this  purpose.  The  first 
model  type,  generally  referred  to  in  the  literature  as  the  pseudo- 
homogeneous  film  model  [5-9],  assumes  that  the  catalyst  layer 
is  a  porous  matrix  comprised  of  Nation,  platinum,  and  carbon  in 
random  (homogeneous)  configuration.  This  model  allows  for  path¬ 
ways  of  gases,  electrons,  and  protons  within  the  catalyst  layer,  and 
captures  some  of  the  essential  transport  phenomena  prevalent  in 
the  catalyst  layer.  However,  this  model  does  not  acknowledge  the 
necessity  for  the  existence  of  the  triple-phase  boundary  for  a  func¬ 
tioning  catalyst  layer.  In  contrast,  the  flooded  agglomerate  concept, 
proposed  in  the  late  1980s  [10,11],  contends  that  the  platinum  is 
supported  on  carbon  particles,  which  forms  agglomerates  when 
mixed  with  an  ionomer.  The  agglomerate  may  even  be  coated  fully 
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or  partially  by  an  additional  ionomer  layer.  The  oxygen  finds  its  way 
to  the  platinum  by  first  dissolving  in  the  ionomer,  and  is  consumed 
as  it  transports  to  the  core  of  the  carbon-platinum  aggregate.  From 
the  preceding  discussion,  it  is  clear  that  this  model  guarantees  the 
existence  of  triple-phase  boundaries  as  long  as  sufficient  amounts 
of  the  ionomer  are  present.  Microscopic  observations  [13,20]  sup¬ 
port  the  structure  proposed  by  this  model,  and  recent  calculations 
performed  using  this  model  appear  to  match  experimental  data 
better  than  the  pseudo-homogeneous  film  model  [8,19].  It  is  worth 
noting  that  the  popularity  of  the  flooded  agglomerate  model  is  in 
part  due  to  the  fact  that  modeling  studies  using  this  model  have 
also  shown  that  the  choice  of  spherical  or  cylindrical  agglomerates 
did  not  significantly  alter  the  overall  performance  of  the  PEMFC 
[27],  indicating  that  the  exact  geometrical  shape  of  the  agglomerate 
is,  perhaps,  not  of  great  importance.  This  finding  presents  another 
reason  not  to  pursue  tedious  DNS  studies. 

Most  notable  among  recent  studies  that  have  used  the  flooded 
agglomerate  model  is  the  work  of  Jaouen  et  al.  [15,16,27],  Sun  et 
al.  [18],  and  Secanell  et  al.  [19].  Jaouen  et  al.  have  successfully 
used  this  model  to  predict  the  performance  of  the  cathode  as  a 
function  of  operating  conditions  and  cathode  layer  thicknesses. 
Their  calculations  predicted  experimentally  observed  double  Tafel 
slopes  [16],  attributed  to  local  (within  agglomerate)  mass  trans¬ 
port  limitations.  Flowever,  the  model  was  not  exercised  to  study 
effect  of  other  compositional  variables  such  as  Nation  loading  and 
platinum  loading.  Such  studies  were  performed  by  Sun  et  al.  [18], 
who  predicted  an  optimum  Nation  loading— consistent  with  exper¬ 
imentally  observed  behavior  [3,4].  Studies  performed  by  both  Sun 
et  al.  [18]  and  Secanell  et  al.  [19]  suggests  that  the  performance  of 
the  cathode  improves  monotonically  as  the  agglomerate  radius  is 
decreased.  However,  it  is  believed  [27]  that  excessive  reduction  in 
the  agglomerate  radius  can  lead  to  severe  reduction  in  the  size  of 
the  macropores  and  high  tortuosity.  This,  in  turn,  can  result  in  dete¬ 
rioration  of  the  performance,  especially  at  high  current  densities, 
due  to  global  mass  transport  limitations. 

In  this  study,  we  generalize  the  flooded  agglomerate  model 
by  carefully  assimilating  findings  of  past  studies,  and  extend¬ 
ing  the  core  model  to  include  the  following  effects:  (1)  effect  of 
Nation  loading,  (2)  effect  of  catalyst  (platinum)  loading,  (3)  effect 
of  carbon-platinum  ratio,  (4)  agglomerate  radius,  and  (5)  cathode 
catalyst  layer  thickness.  In  an  effort  to  do  so,  care  has  been  exercised 
to  account  for  both  local  (agglomerate-scale)  as  well  as  global  (cat¬ 
alyst  layer  thickness-scale)  effects.  Whenever  possible,  an  attempt 
has  also  been  made  to  compare  predicted  results  with  experimental 
findings,  at  least  qualitatively. 


2.  Model  description 


Fig.  1.  Schematic  representation  of  a  cross-section  of  a  PEM  fuel  cell.  The  magni¬ 
fied  portion  of  the  cathode  catalyst  layer  shows  a  representation  of  the  uniformly 
distributed  spherical  agglomerates  surrounded  by  gas  pores. 


solution  into  the  polymer  coating  followed  by  diffusion  through 
the  polymer  within  the  agglomerate.  The  purpose  of  the  sub-grid 
scale  agglomerate  model  is  to  provide  a  means  to  calculate  the 
volumetric  current  density  (Am-3)  in  the  cathode  catalyst  layer. 
Expressions  to  this  effect  are  derived  in  this  section.  Expressions 
for  volume  fractions  of  the  three  components  of  the  cathode,  which 
influence  various  effective  transport  properties,  are  also  derived  in 
this  section. 

The  oxygen  reduction  reaction  (ORR)  kinetics  in  the  polymer 
electrolyte  is  assumed  to  follow  a  Tafel  law  and  to  be  first-order 
in  oxygen  concentration.  Hence,  the  surface  current  density  on  a 
smooth  catalyst  surface  is  given  by  [15]: 


where  is  the  standard  reference  oxygen  concentration,  ij*f  is  the 
standard  exchange  current  density  for  a  smooth  catalyst  surface, 
and  r\  is  the  local  overpotential  defined  as  the  difference  between 
the  electronic  and  protonic  phase  potentials,  i.e.,  r]  =  </>s  -  </>f-  c  and  c* 
are  the  dissolved  oxygen  concentration  in  the  polymer  at  the  cata¬ 
lyst  surface  and  the  dissolved  oxygen  concentration  in  the  polymer 
in  equilibrium  with  the  inlet  gas,  respectively. 

Within  the  agglomerate,  oxygen  is  assumed  to  diffuse  through 
the  polymer  electrolyte  to  the  catalyst  surface  according  to  Fields 
law.  While  it  diffuses  through  the  agglomerate,  it  reacts  on  the 
platinum  clusters  impregnated  on  the  carbon  surfaces,  and  is  con¬ 
sumed  (Fig.  2).  This  phenomenon  is  best  described  by  the  following 
reaction-diffusion  balance  equation  [15]: 


D 


eff 

02,N 


/ d2c  2  dc\ 
y  dr2  +  r  dr  J 


A4'sat(l  ~  £agg) 

r  nF 


=  0 


y  =  1  for  0  <  r  <  ragg;  y  =  0  for  ragg  <  r  <  (ragg  +  5)  (2) 


2  A.  Agglomerate  model 


The  boundary  conditions  for  the  above  equation  are  as  follows  [15]: 


In  our  agglomerate  model,  the  cathode  is  assumed  to  be  com¬ 
prised  of  a  large  number  of  spherical  agglomerates  surrounded  by 
gas  pores  (Fig.  1 ).  An  agglomerate,  in  turn,  is  comprised  of  clusters 
of  carbon  supporting  platinum  particles  held  together  by  proton¬ 
conducting  polymer.  The  agglomerate  may  be  coated  with  a  thin 
film  of  polymer  (Fig.  2).  Thus,  the  cathode  consists  of  three  compo¬ 
nents:  the  solid  carbon-platinum  clusters,  the  polymer  electrolyte, 
and  gas  pores.  For  the  model  to  be  valid,  the  agglomerate  size  must 
be  sufficiently  small  compared  to  the  cathode  layer  thickness  such 
that  the  potential  within  an  agglomerate  may  be  assumed  to  be 
constant.  Transport  of  reactants  from  the  cathode  channel  to  the 
cathode  active  layer  is  dominated  by  diffusion  within  the  pores 
of  the  cathode  diffusion  layer  and  the  cathode  catalyst  layer  (see 
Fig.  1 ).  Oxygen  reaches  platinum  within  the  agglomerate  by  dis- 


dc 

dr 


=  0 


(3a) 


Agglomerate  Nucleus 


Polymer  Coating 


Carbon-supported 

Catalyst 


Fig.  2.  Agglomerate  nucleus  with  carbon-supported  catalyst  particles  distributed  in 
polymer  with  a  polymer  coating  around  the  nucleus. 
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c(ragg)  —  c(ragg) 


'agg 

neff  dc : 
^02, N  dr 


dc 

=  D°2,N3? 


c(ragg  +  5)  -  —  Co2,g  -  Hco2,g 
co2,g 


(3b) 

(3c) 

(3d) 


where  ragg  is  the  radius  of  the  agglomerate  nucleus,  8  is  the 
thickness  of  the  polymer  coating  around  the  agglomerate,  £agg  is 
the  volume  fraction  of  polymer  present  within  the  agglomerate 
nucleus,  andAv  is  the  total  catalyst  surface  area  per  unit  volume.  Eq. 
(3d)  represents  equilibrium  between  dissolved  oxygen  concentra¬ 
tion  in  Nation  (c)  and  oxygen  concentration  in  the  gas  pores  (co2,g) 
at  the  outer  surface  of  the  polymer  coating,  where  H  is  Henry’s 
constant.  The  term  c*2  g  in  Eq.  (3d)  represents  the  concentration 

of  oxygen  at  cathode  inlet.  N  is  the  effective  oxygen  diffusion 
coefficient  in  Nation  within  the  agglomerate,  and  is  obtained  from 
the  Bruggemann  relation  as  [28]: 

d02,n  =  Do2,Nfilg|  for  r  <  ragg 
<N  =  do2,n  f°r  ragg  <  r  <  (ragg  +  5) 


The  system  of  equations,  Eqs.  (l)-(3),  can  be  solved  analytically. 
It  yields  an  expression  for  the  volumetric  current  density  in  the 
cathode,  written  as  [15]: 

(\  l-ocj/n 


x  exp 


Effi  Eff2(l  -  £Cat) 


agg 


Co 


'2>g 


(ragg  +  <5)3  co2,g 


(5) 


where  £cat  is  the  porosity  of  the  cathode  catalyst  layer,  c02,g  and 
cj2  g  are  the  oxygen  concentrations  at  the  polymer-coating-gas- 
pore  interface  and  at  the  inlet,  respectively.  Ef^  and  Eff2  are  the 
nucleus  and  film  effectiveness  factors.  They  represent  corrections 
to  the  pure  kinetic  current  density  due  to  diffusion  limitations  in 
the  agglomerate  nucleus  and  the  polymer  coating,  respectively.  Eff[ 
and  Eff2  are  written  as  follows  [15]: 


Effi  = 


3 

9ragg 


(  1 

V  tanh(  qragg) 


1 

Qrzgg 


(6) 


1  +  (l/3)(ra2gg5)/(ragg  +  Stfel ^Eff, 

where 

2  Av(l  -  eagg) (c/cg)  W"  exp  (-faF)/ (RT)t,) 
nFDo2,Nc*£agg 


Since  Effi  and  Eff2  represent  corrections  to  the  pure  kinetic  cur¬ 
rent  density,  it  is  informative  to  investigate  how  they  change  with 
ragg,  sagg  and  8  in  order  to  gain  an  understanding  of  how  the 
agglomerate-level  mass  transport  losses  behave.  Fig.  3  shows  a 
plot  of  these  two  quantities.  It  is  seen  from  Fig.  3(a)-(c)  that  the 
local  mass  transport  losses  are  largest  (Eff[  is  smallest)  when  the 
agglomerate  radius,  ragg,  is  large  and  the  Nation  volume  fraction 
within  the  agglomerate,  £agg,  is  small.  And  the  losses  are  larger  at 
larger  values  of  the  cathode  overpotential.  It  is  also  observed  that 
the  film  coating  has  almost  no  effect  on  local  mass-transport  losses, 
probably  because  the  film  coating  thickness  is  small  in  the  first 
place. 

The  parameter  Av  represents  the  total  catalyst  surface  area  per 
unit  volume  of  cathode  catalyst  layer  available  for  the  ORR.  This 
value  is  related  to  the  platinum  mass  loading  (m pt),  the  platinum 


particle  size  and  the  platinum  to  carbon  mass  ratio  (Pt|C)  in  the  cat¬ 
alyst  layer  ink.  The  dependence  of  this  parameter  on  the  platinum 
mass  loading  is  given  as  [19]: 


AV=A0^ 


(9) 


where  A0  is  the  total  catalyst  surface  area  per  unit  mass  of  catalyst 
particle.  Based  on  empirical  data  provided  by  Marr  and  Li  [6]  for 
various  platinum  supported  catalysts,  Secanell  et  al.  [19]  estimated 
this  value  using  the  following  correlation: 

Ao  =  2.2779  x  106(Pt|C)3  -  1.5857  x  106(Pt|C)2 

-2.0153  x  106(Pt|C)  +  1.5950  x  106  (10) 


In  order  to  compute  the  effective  properties  for  transport  in  the 
cathode  catalyst  layer,  the  volume  fractions  of  the  three  compo¬ 
nents  of  the  cathode,  namely  the  solid  component  comprised  of 
carbon-platinum  clusters,  the  proton-conducting  polymer  compo¬ 
nent  and  the  porous  component  for  mass  transport,  need  to  be 
calculated.  The  composition  of  the  catalyst  layer  ink  used  strongly 
influences  these  volume  fractions. 

The  solid  component  volume  fraction  is  given  by  the  amounts  of 
platinum  and  carbon  in  the  cathode  layer.  These  values  are  obtained 
from  the  platinum  mass  loading  (mPt)  and  the  platinum  to  car¬ 
bon  mass  ratio  (Pt|C)  in  the  catalyst  layer  ink.  The  solid  component 
volume  fraction  is  obtained  as  [19]: 


(J_  1  -Pt|C\  rn pt 
vPpt  Pt|Cpc  )  L 


(ID 


where  pPt  and  pc  are  the  densities  of  platinum  and  carbon  respec¬ 
tively,  and  L  is  the  cathode  catalyst  layer  thickness. 

Under  the  assumption  that  the  cathode  catalyst  layer  consists 
of  spherical  agglomerates,  and  that  the  agglomerates  are  made 
up  of  only  the  solid  component  and  the  polymer,  the  number  of 
agglomerates  per  unit  volume  required  to  obtain  the  above  solid 
component  volume  fraction  is  given  by  the  following  equation 
[19]: 


(4/3)7rragg(l  —  £agg) 


(12) 


Assuming  that  the  entire  polymer  electrolyte  in  the  cathode  catalyst 
layer  is  only  present  either  in  the  agglomerate  nucleus  or  in  the 
polymer  coating  around  the  agglomerate,  the  volume  fraction  of 
the  polymer  in  the  cathode  can  be  obtained  from  [19]: 


£n  -  o  7th  [rfgg£agg  +  {(ragg  +  <$)3  -  rfgg  ] 


(13) 


where  8  is  the  average  thickness  of  the  polymer  film  coating  around 
the  agglomerate. 

Once  the  volume  fractions  of  the  solid  component  and  the  poly¬ 
mer  are  obtained,  the  porosity  of  the  cathode  catalyst  layer  is 
obtained  as 


£cat  =  l-£s-£N  (14) 

This  volume  fraction  represents  the  volume  of  macro-pores  in 
the  cathode  catalyst  layer  which  is  available  for  transport  of 
reactants. 

As  can  be  seen  from  Eqs.  (1),  (2),  (3a)-(3d),  (4)-(14),  the  trans¬ 
port  properties,  as  well  as  the  parameters  for  electrochemical 
reaction  kinetics  are  simultaneously  influenced  by  the  radius  of 
agglomerate  nucleus  (ragg),  the  agglomerate  coating  thickness  (5), 
volume  fraction  of  polymer  within  the  agglomerate  nucleus  (£agg), 
platinum  mass  loading  (mPt),  platinum-carbon  mass  ratioin  the 
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Fig.  3.  Nucleus  effectiveness  factor  (Effi )  with  varying  agglomerate  radius  and  polymer  volumetric  fraction  within  agglomerate  nucleus  for  nominal  cathode  overpotential 
values  of  (a)  r]c at  =  -0.3  V,  (b)  ijc at  =  -0.4  V,  and  (c)  ijc at  =  -0.5  V.  Using  the  maximum  value  of  the  nucleus  effectiveness  factor,  the  effect  of  the  polymer  coating  thickness  (5) 
on  the  film  effectiveness  factor  (Eff2)  is  presented  in  (d).  Kinetic  constants  for  these  calculations  were  obtained  using  L  =  15  |jim,  mPt  =  0.4  mg  cm-2  and  Pt|C  =  0.28. 


catalyst  layer  ink  (Pt|C),  and  the  cathode  catalyst  layer  thickness 

(L). 

2.2.  Computational  fluid  dynamics  (CFD)  model 

The  sub-grid  scale  catalyst  layer  model,  described  in  the 
preceding  section,  was  implemented  into  a  two-dimensional  com¬ 
putational  fluid  dynamics  (CFD)  code  that  has  been  developed 
specifically  for  PEMFC  applications,  and  has  been  validated  [29]. 
The  purpose  of  this  step  is  to  enable  prediction  of  the  overall  perfor¬ 
mance  (polarization  behavior)  of  the  fuel  cell  by  taking  into  account 
both  local  and  global  losses.  The  governing  equations  are  the  equa¬ 
tions  of  conservation  of  mass  (both  overall  and  individual  species), 
momentum,  energy,  and  charge.  For  simplicity,  it  is  assumed  that 
the  temperature  does  not  change  spatially,  and  thus,  the  energy 
conservation  equation  is  not  solved.  Furthermore,  it  is  assumed 
that  water  exists  only  in  its  vapor  phase  in  gas  diffusion  layers 
and  channels,  and  two-phase  effects  are  not  considered.  The  gov¬ 
erning  conservation  equations  are  different  inside  and  outside  the 
membrane  and  are,  therefore,  presented  here  separately. 

2.2.1.  Channels ,  gas  diffusion  layers  (GDL),  and  active  catalyst 
layers  (ACL) 

The  governing  conservation  equations  for  mass  and  momentum, 
are  written  as  [30-32]: 

9 

overall  mass  :  — ( sp )  +  V  •  (epU)  =  5m  (15) 

at 


momentum  :  ^-(peU)+V  •  (peUU)=  -  eVp+V  •  (/teffVU)  + 

at  K 

(16) 

f) 

species  mass  :  ^(pYk)+V  •  {p\JYk)=  -  V  J k+Sk  Vk  =  1, 2,  . . . ,  N 

(17) 

where  p  is  the  density,  p  is  the  pressure,  /z  is  the  dynamic  viscosity 
of  the  fluid,  and  U  is  the  fluid  velocity  vector,  s  is  the  porosity  of 
the  medium  and  /c  is  the  permeability.  In  Eq.  (15),  5m  is  the  mass 
source  term  due  to  electrochemical  reactions,  and  is  non-zero  only 
in  the  catalyst  layers.  It  results  from  the  fact  that  a  conservation 
equation  for  protons  is  not  directly  solved,  and  therefore,  the  mass 
of  protons  created  or  destroyed  need  to  be  subtracted  out  of  the 
overall  continuity  equation.  The  last  term  in  Eq.  (16)  represents  the 
sub-grid  scale  drag  force  imposed  by  the  pore  walls  on  the  fluid, 
written  using  the  linear  Darcy’s  Law  [28].  In  purely  open  regions, 
such  as  in  the  gas  channels,  e  ->  1  and  /c  ->►  oo,  and  Eqs.  (15)  and  (16) 
reduce  to  the  well-known  Navier-Stokes  equations.  In  Eq.  (17),  Yk 
is  the  mass  fraction  of  the  kth  species,  Jk  is  the  mass  diffusion  flux 
of  the  kth  species,  and  Sk  is  the  production  rate  of  the  kth  species 
due  to  electrochemical  reactions.  The  total  number  of  species  in  the 
system  is  denoted  by  N. 

The  diffusion  flux  of  the  kth  species,  Jk ,  is  modeled  using  the 
so-called  dilute  approximation  [33].  By  this  approximation,  the 
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diffusion  flux  is  written  as 

ik  =  -pDkmVYk  (18) 

where  Dkm  is  the  free-stream  diffusivity  of  species  k  into  the  mix¬ 
ture,  and  is  henceforth  denoted  by  Dk  for  simplicity.  The  free-stream 
diffusivity  is  given  by  the  following  relation  [33]: 

Dkm  =  Dk=  l~X\  (19) 

Z^i=l,i/k  Dki 

where  Dki  is  the  binary  diffusivity  of  species  k  into  species  i,  and  Xk 
is  the  mole  fraction  of  the  kth  species.  Substitution  of  Eq.  (18)  into 
Eq.  ( 17)  yields  the  appropriate  species  transport  equation  under  the 
dilute  approximation  formulation  for  open  regions: 

| XspYk)+V-(epUYk)=V-(pDkVYk)+Sk  Vk  =  1,  2, . . . ,  N  (20) 

In  porous  regions,  it  is  customary  to  use  the  so-called  Bruggemann 
relation  [28],  and  modify  the  free-stream  diffusion  coefficient,  such 
that  the  governing  species  conservation  equation  becomes 

^(epVk)  +  V  .  (epUYfc)  =  V  .  ( pDkeWYk )  +  Sk  Vk  =  1 , 2, . . . ,  N 

(21) 


Fig.  4.  Effect  of  agglomerate  radius  (ragg)  on  non-dimensional  effective  gas  dif- 
fusivities  in  cathode  for  L  =  15|jim,  <5/ragg  =  0.06,  £agg  =  0.4,  mPt  =  0.45  mg  cm-2  and 
Pt|C  =  0.28.  In  the  figure,  D^ff  stands  for  the  effective  diffusivity  of  species  k  into 
the  gas  mixture  in  the  cathode  catalyst  layer,  while  D£at  stands  for  the  free  stream 
diffusivity  of  species  k  into  the  gas  mixture  in  the  cathode  catalyst  layer. 


where  r  is  the  tortuosity  of  the  porous  region.  Traditionally,  a  value 
of  r  =  15  is  used  in  Eq.  (21 ).  In  order  to  examine  the  effect  of  cathode 
structure  and  tortuosity  on  fuel  cell  performance,  a  more  sophis¬ 
ticated  tortuosity  model,  proposed  by  Abbasi  et  al.  [34],  is  used. 
Hence  the  governing  species  conservation  equation  in  the  cathode 
catalyst  layer  becomes: 

|(ficat  pYk)  +  V  .  (ecatpUYfc)  =  V  .  (pDfcgvYfc) 

+  Sk  Vk  =  1, 2, . . . ,  JV  (22) 


where  the  tortousity  rcat  is  given  by  [34]: 

teat  =  —  +  1.196^  (23) 

£cat  d 

where  <rdev  is  the  standard  deviation  of  the  pore  size  in  the  cathode 
catalyst  layer.  In  this  study,  a  constant  value  of  crdev  =  50  nm  has 
been  used  for  the  sake  of  convenience,  d  is  the  average  pore  size  of 
the  cathode  catalyst  layer,  and  is  given  by  [35]: 


4  gcat 

3(1—  £cat) 


(24) 


Clearly,  this  tortuosity  model  explicitly  accounts  for  the  agglomer¬ 
ate  radius  variation  in  addition  to  variations  in  the  global  porosity 
(as  stipulated  by  the  Bruggemann  model).  Thus,  one  can  account 
for  changes  in  tortuosity  (and  mass  transport  resistance)  due  to 
change  in  the  agglomerate  size  while  keeping  the  overall  porosity 
unchanged— an  effect  that  cannot  be  predicted  by  the  Bruggemann 
model.  The  differences  between  the  diffusion  coefficients  predicted 
by  the  two  models  are  shown  in  Fig.  4.  The  key  observation  from  this 
figure  is  that  for  a  constant  porosity,  the  Bruggemann  model  does 
not  predict  any  variation  of  the  diffusion  coefficient  with  agglomer¬ 
ate  (or  pore)  size,  while  the  model  used  here  predicts  a  sharp  drop 
in  the  effective  diffusion  coefficient  when  the  agglomerate  (or  pore) 
size  becomes  significantly  small. 

The  source  due  to  the  electrochemical  reactions  is  non-zero  only 
in  the  active  catalyst  layers  of  the  anode  and  cathode,  and  is  zero 
elsewhere.  It  is  written  as  [36]: 


Sk  = 


jjMk 
2  nF 


(25) 


of  electrons  transferred  during  the  electrochemical  reaction,  and  F 
is  the  Faraday  constant.  (S/V)eff  represents  the  ratio  of  the  active 
surface  area  of  the  catalyst  to  its  overall  volume  (including  the 
support).  The  anode  transfer  current  jf1  is  expressed  using  the 
Butler-Volmer  equation  [36]: 


J?  =J  o 


exp(w")-‘ 


“p  (-£")]  ni'W 

k= 1 


Pk 


(26) 


where  j0  is  the  reference  current  density,  aa  and  ac  are  kinetic 
constants  determined  from  experimentally  measured  Tafel  plots. 
[Ak]  and  jdk  are  the  molar  concentrations  and  the  concentra¬ 
tion  exponents  for  the  kth  species,  respectively,  rj  is  the  electrode 
overpotential  (including  both  activation  and  concentration  overpo¬ 
tential)  and  is  defined  as  the  difference  between  the  electronic  or 
solid  phase  potential  (</>s)  and  the  electrolyte  or  ionic  phase  poten¬ 
tial  (0F),  i.e.,  t)  =  0s  —  0f*  In  this  study,  the  above  Butler-Volmer 
kinetic  relation  was  used  only  for  the  anode,  while  the  cathode 
transfer  current  was  obtained  from  Eq.  (5)  (Section  2.1 ). 

In  a  PEMFC  calculation,  in  addition  to  mass  conservation,  it  is 
necessary  to  enforce  charge  conservation.  Under  the  assumption  of 
electro-neutrality,  charge  conservation  reduces  to  current  conser¬ 
vation,  written  as 


V  •  i  =  0  (27) 

where  i  is  the  current  density  vector.  In  a  PEMFC,  the  current  flow  is 
due  to  protons  (H+)  flowing  through  the  membrane,  resulting  in  an 
ionic  phase  current  (iF),  and  due  to  electrons  flowing  through  the 
carbon  in  the  porous  matrix  of  the  gas  diffusion  layers,  resulting  in 
an  electronic  phase  current  (is).  Thus,  Eq.  (27)  can  be  rewritten  as 
follows: 


V  •  iF  +  V  •  is  =  0  (28) 

Since  current  transport  in  the  ionic  phase  is  due  to  ions  and  that  in 
the  electronic  phase  is  due  to  electrons,  the  transport  in  each  phase 
is  governed  by  separate  electric  potential  fields.  Using  Ohm’s  law, 
Eq.  (28)  can  be  written  as 

V  •  (<7FV0F)  +  V  •  (<tsV0s)  =  0  (29) 


where  jj  is  the  net  transfer  current  due  to  electrochemical  reac-  where  <rF  and  o $  are  the  conductivities  of  the  ionic  and  electronic 
tion,  Mk  is  the  molecular  weight  of  the  kth  species,  n  is  the  number  phases,  respectively.  The  exchange  of  current  from  the  ionic  to  the 
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electronic  phase  occurs  due  to  electrochemical  reactions  during 
which  electrons  are  transferred  from  one  phase  to  the  other.  Thus, 
Eq.  (29)  can  be  rewritten  as  [31,36]: 

-V  .  (ofV</>f)  =  V  .  (orsV</>s)  =  >r  (£)  (30) 

Eq.  (30)  represents  two  elliptic  partial  differential  equations  that 
are  strongly  coupled  through  the  transfer  current  source.  The  ionic 
phase  electric  potential  equation  (for  0F)  must  be  solved  in  the  ACL 
and  the  membrane,  while  the  electronic  phase  electric  potential 
equation  (for  0S)  must  be  solved  in  the  ACL  and  the  GDL  (see  Fig.  1 ). 
In  the  ACL,  both  equations  are  solved,  and  are  strongly  coupled.  The 
difference  in  value  between  0S  and  4>p  represents  the  total  electrode 
overpotential. 


2.2.2.  Water  and  current  transport  in  membrane 

Nation  membranes  are,  for  all  practical  purposes,  impermeable 
to  all  gases  except  water.  Thus,  the  governing  equation  for  species 
transport  outside  the  membrane  (Eq.  (17))  is  irrelevant  in  this  con¬ 
text.  The  only  species  that  needs  to  be  considered  is  water.  For  all 
other  species,  zero  flux  boundary  conditions  must  be  used  at  the 
membrane-ACL  interface.  Water  transport  in  the  membrane  of  a 
PEMFC  occurs  primarily  due  to  diffusion  and  electro-osmotic  drag, 
although  there  is  some  evidence  that  pressure-driven  advection 
may  also  occur  [37-41].  In  phenomenological  membrane  models, 
it  is  customary  to  express  water  transport  in  terms  of  the  water 
content,  A,  as  [42]: 


V° 


=  v° 


Pm 

Mm 


D>VA 


(31) 


The  water  content  is  defined  as  the  ratio  of  the  number  of  water 
molecules  to  the  number  of  S03_H+  charge  sites  in  the  medium.  In 
Eq.  (31 ),  it  has  been  assumed  that  advective  transport  is  negligible. 
Dx  denotes  the  diffusion  coefficient  of  water  in  Nation  expressed 
in  terms  of  the  water  content,  pd  denotes  the  electro-osmotic  drag 
coefficient,  and  pm/Mm  represents  the  molar  density  of  the  mem¬ 
brane.  Fuller  and  Newman  [43]  give  the  diffusion  coefficient  as: 


Fig.  5.  Effect  of  polymer  coating  thickness  (5)  on  non-dimensional  effective  protonic 
conductivity  of  cathode  for  L  =  15fxm,  ragg  =  1000nm,  £agg  =  0.4,  mPt  =  0.4mgcirr2 
and  Pt|C  =  0.28. 


ductivity  in  the  cathode  is  given  by  Jaouen  et  al.  [15]  as: 


°rpat  —  (1  —  £cat) 


(gagg  ~  l) 

(l  +  <Vragg) 


tTp 


(36) 


Proton  transport  between  agglomerates  requires  that  there  be 
sufficient  polymer  coating  around  agglomerates.  As  the  coating 
thickness  decreases,  the  contact  between  agglomerates  decreases, 
which  in  turn  leads  to  decreased  proton  transport.  To  better  capture 
this  trend,  we  need  that  as  8  ->  0,  <7pat  -*  0  as  well.  To  this  end,  in 
the  present  model,  the  effective  ionic  phase  electrical  conductivity 
in  the  cathode  is  modified  to 


—  (1  —  £cat) 


1  +  ■ 


(gagg  -  1 ) 


(l  +  <$/ragg  +  Clo) 


Op 


(37) 


Dx  =  (2.1  x  10-7)  x  A  x  exp - — 

The  electro-osmotic  drag  coefficient  is  given  by  [44]: 


(32) 


where  a0  is  given  by: 


a0  =  min 


0, 


(38) 


ru  =  (33) 

Current  transport  in  the  membrane  is  described  by  the  ionic 
phase  part  of  Eq.  (30),  except  that  there  is  no  source,  i.e.,  the  right 
hand  side  of  the  equation  is  zero.  The  ionic  phase  electrical  con¬ 
ductivity  of  the  membrane,  op,  is  again  expressed  by  the  empirical 
correlation  [44]: 


op  =  100  exp 


,268  (s® 


(34) 


where  030  is  the  electrical  conductivity  of  the  membrane  at  30  °C, 
and  is  given  by  [43]: 


ct30  =  0.005139A.  -  0.00326  for  X  >  1 


(35) 


The  ionic  phase  electrical  conductivities  in  the  anode  catalyst 
layer  and  the  cathode  catalyst  layer  are  evaluated  using  Eqs.  (34) 
and  (35)  based  on  local  values  of  water  content.  However,  in  these 
regions,  the  conductivities  need  to  be  corrected  in  order  to  account 
for  the  fact  that  the  proton-conducting  polymer  is  just  one  of  the 
components  in  the  layer.  The  effective  ionic  phase  electrical  con- 


Eq.  (37)  satisfies  the  requirement  that  very  low  values  of  8  lead  to 
very  low  values  of  <7pat,  while  returning  the  values  obtained  from 
Eq.  (36)  at  sufficiently  high  values  of  8,  as  shown  in  Fig.  5.  It  is 
worth  noting  that  the  contention  that  the  agglomerate  is  coated 
with  a  uniformly  thick  layer  of  the  ionomer  is  an  assumption  in  the 
flooded  agglomerate  model.  Clearly,  that  would  make  it  impossible 
for  electrons  to  flow  from  one  carbon  particle  to  another.  However, 
in  a  real  3D  microstructure,  there  is  always  some  partial  contact 
between  the  cores  of  adjacent  agglomerates,  which  makes  electron 
transport  possible. 

It  is  evident  from  Eqs.  (31)-(35)  that  the  two  quantities  that 
dictate  the  mass  transport  and  electrical  properties  of  a  Nation 
membrane  are  its  water  content  and  temperature.  The  water  con¬ 
tent  is  obtained  by  solution  of  the  conservation  equation  for  water 
(Eq.  (31)).  However,  solution  of  this  equation  requires  specifica¬ 
tion  of  either  the  water  content  itself  or  the  flux  of  water  at  the 
membrane-ACL  interface,  which  in  turn  requires  coupling  with 
the  rest  of  the  computational  domain.  Various  approaches  for  cou¬ 
pling  the  membrane  with  the  overall  calculation  procedure  and 
for  implementation  of  the  interface  conditions  are  discussed  in 
Kamarajugadda  and  Mazumder  [29]. 
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3.  Results  and  discussion 

The  governing  conservation  equations  of  mass,  momentum  and 
current,  described  in  the  preceding  section,  were  solved  using 
a  conservative  finite-volume  technique  [45].  The  SIMPLE  algo¬ 
rithm  [45]  was  used  to  address  pressure-velocity  coupling  in  the 
Navier-Stokes  equations.  The  2D  model,  shown  in  Fig.  1,  was  used 
for  simulations.  All  simulations  were  performed  on  a  uniform  grid 
with  50  cells  in  the  axial  direction  and  150  cells  in  the  cross-flow  (y) 
direction.  Nominally,  40  cells  were  used  across  each  flow  channel, 
10  across  each  GDL,  5  across  each  active  catalyst  layer,  and  40  across 
the  membrane.  This  particular  mesh  size  was  chosen  after  a  rigor¬ 
ous  grid-independence  study,  the  details  of  which  can  be  found  in 
Kamarajugadda  and  Mazumder  [29]. 

Moist  hydrogen  and  moist  air  was  introduced  into  the  anode 
and  cathode  inlet,  respectively.  A  plug  velocity  profile  was  imposed, 
and  the  magnitude  of  the  velocity  was  calculated  from  a  prescribed 
equivalent  current  density.  For  all  simulations,  the  temperature  was 
assumed  to  be  50  °C  everywhere,  and  the  energy  equation  was  not 
solved.  Other  relevant  parameters  used  for  simulations  are  reported 
in  Table  1 .  The  kinetic  constant,  j0,  for  the  electrochemical  reactions 
at  the  anode  was  determined  by  calibrating  the  current  density  at 
low  bias  voltage  against  experimental  data  reported  by  Ticianelli  et 
al.  [46,47]— a  practice  that  has  been  used  in  the  past  [31  ].  All  simu¬ 
lations  were  run  with  a  prescribed  bias  voltage  boundary  condition 
rather  than  a  prescribed  current  boundary  condition. 

All  transport  properties  of  the  fluid,  namely  viscosity  and  binary 
diffusion  coefficients  were  computed  using  the  Chapman-Enskog 
equations  of  kinetic  theory  [33,48].  The  Lennard-Jones  potentials, 
which  are  needed  as  inputs,  were  obtained  from  the  CHEMKIN 
database.  The  density  of  the  fluid  was  calculated  using  the  ideal  gas 
law.  The  solutions  were  deemed  to  be  converged  when  the  residuals 
of  each  of  the  equations  decreased  by  seven  orders  of  magnitude. 
In  order  to  compute  the  actual  cell  voltage  from  the  prescribed  bias 
voltage  (or  potential  loss),  it  is  necessary  to  know  the  open  circuit 
potential.  As  done  by  previous  researchers  [31,49],  rather  than  use 


Table  1 

Simulation  parameters  and  values  of  key  properties  at  323 1<  and  1/1  atm 


Model  parameter 

Value/method  of  calculation 

Gas  channel  length 

0.07112  m 

Gas  channel  width 

7.62  x  10-4  m 

Diffusion  layer  width 

2.54  x  10“4  m 

Membrane  width 

1.75  x  10-4  m 

Membrane  permeability 

1.8  x  10“18  m2 

Diffuser  and  catalyst  layer 

1.76  x  10-11  m2 

permeability 

Membrane  porosity 

0.28 

Anode  and  cathode  diffuser 

0.5 

layer  porosity 

Tortuosity  for  Bruggemann 

1.5 

correlation 

Air  side/fuel  side  pressures 

1/1  atm 

Relative  humidity  of  inlet 

100% 

streams 

Air  side  inlet  N2/O2  molar 

79/21 

ratio 

H2  stoichiometric  flow 

2.8  A  cm-2  equiv. 

O2  stoichiometric  flow 

3.0  A  cm-2  equiv. 

Transfer  coefficients  at 

0.5 

anode  (Tafel  constants) 

Concentration  dependence 

0.5  (H2) 

at  anode 

Reference  current  density  at 

1.6  x  1011  (A  nrr3 )  (m3  kmol  H2 )_1/2 

anode 

Membrane  electrical 

Springer  et  al.  [1] 

conductivity 

Diffuser  layer  electrical 

100  (ISm-1) 

conductivity 

Table  2 

Values  of  cathode  model  parameters  and  properties  for  baseline  case 


Model  parameter 

Value/method  of  calculation 

Cathode  Layer  thickness 

15  [Jim 

Agglomerate  radius 

1000  nm 

Polymer  coating  thickness 

80  nm 

Polymer  volumetric  fraction  in 

0.4 

agglomerate  nucleus 

Platinum  mass  loading 

0.4  mg  cm  2 

Platinum-carbon  ratio  in  catalyst 

0.28 

layer  ink 

Transfer  coefficient  at  cathode 

1  (>0.8  V)  [18] 

(Tafel  constant) 

Reference  exchange  current 

0.55  (<0.8  V)  [18] 

1.288  x  10“4  Acirr2  (>0.8  V)  [18] 

density  at  cathode 

Diffusion  coefficient  of  02  in  Nation 

2.291  x  10-2  Acirr2  (<0.8  V)  [18] 
3.966  x  10-10  m2  s-1  [18] 

Henry’s  constant 

10  [18] 

the  Nernst  potential,  an  empirical  correlation  [50]  was  used  for  the 
open  circuit  potential: 

Voc  =  0.0025T  +  0.2329  (39) 

This  correlation  results  in  an  open  circuit  voltage  of  1.04  V  at  50  °C. 
The  model  was  first  validated  (Section  3.1)  against  results  from 
Jaouen  et  al.  [15,16]  prior  to  exercising  it  to  study  the  effect  of 
the  various  structural  and  compositional  parameters  of  the  cath¬ 
ode  catalyst  layer  on  fuel  cell  performance.  The  effects  of  cathode 
composition  and  cathode  structure  are  presented  separately  in  Sec¬ 
tions  3.2  and  3.3.  Baseline  values  for  cathode  model  parameters  and 
properties  are  listed  in  Table  2. 

3.1.  Model  validation 

The  cathode  model  described  in  Section  2.1  was  incorporated 
into  the  CFD  model,  described  in  Section  2.2,  developed  for  pre¬ 
dicting  the  PEMFC  performance.  A  two-stage  validation  study  was 
performed,  in  which  the  polarization  curves  were  obtained  for  con¬ 
ditions  described  in  Jaouen  et  al.  for  a  theoretical  baseline  case 
[15]  and  a  given  experimental  case  [16].  In  the  first  stage,  the  per¬ 
formance  of  the  cathode  catalyst  model  was  compared  with  the 
baseline  case  described  in  Table  1  of  Ref.  [15]  while  minimizing  the 
mass  transport  losses  in  the  gas  diffusion  layers,  Ohmic  losses  in  the 
membrane,  and  activation  and  mass  transport  losses  in  the  anode. 
In  this  case,  the  model  was  run  for  a  fuel  cell  operating  at  50  °C  with 
H2/O2  gas  streams  at  1  atm  with  70%  relative  humidity  in  the  anode 
and  cathode  channels.  The  cathode  catalyst  layer  was  assumed  to 
be  L  =  10[xm  thick  with  an  agglomerate  radius  of  ragg  =  0.5  f±m,  a 
polymer  coating  thickness  of  8  =  0  |jim,  and  a  polymer  volumetric 
fraction  of  £agg  =  0.3.  The  porosity  of  the  cathode  catalyst  layer  was 
assumed  to  be  £Cat  =  0.3.  Other  details  of  the  model  are  described 
in  Table  1  of  Ref.  [15].  Fig.  6(a)  shows  the  comparison  between 
the  performance  predicted  by  the  present  model  and  the  perfor¬ 
mance  predicted  by  Jaouen  et  al.  The  double  Tafel  slope  predicted 
by  Jaouen  et  al.  [15]  is  reproduced  in  the  present  model,  and  the  dif¬ 
ference  between  the  current  densities  predicted  by  the  two  models 
is  insignificant. 

In  the  second  stage  of  the  validation  test,  the  model  developed 
was  compared  against  experimental  results  provided  in  Fig.  7  of 
Ref.  [16]  for  Cathode  C.  In  this  case,  the  model  was  run  for  a  fuel  cell 
operating  at  50  °C  with  H2/02  gas  streams  at  1  atm  with  86%  rela¬ 
tive  humidity  in  the  anode  and  cathode  channels.  Dimensions  for 
Cathode  C  with  a  thickness  of  1  =  4.7  p,m,  as  described  in  Ref.  [16], 
were  used,  with  an  agglomerate  radius  of  ragg  =  0.11  pum,  a  polymer 
coating  thickness  of  8  =  0  |jim,  and  a  polymer  volumetric  fraction 
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Fig.  6.  Model  Validation  Studies:  (a)  comparison  with  model  baseline  case  from 
Table  1  of  Jaouen  et  al.  [15]  and  (b)  comparison  with  experimental  results  for  “Cath¬ 
ode  C”  from  Fig.  7  of  Jaouen  et  al.  [16]. 


of  £agg  =  0.38.  The  porosity  of  the  cathode  was  £cat  =  0.27.  Other 
details  regarding  the  parameters  used  in  the  model  are  available 
in  Ref.  [16].  Fig.  6(b)  shows  the  comparison  between  the  fuel  cell 
performance  predicted  by  the  current  model  and  the  polarization 
curve  obtained  from  the  experiment.  Once  again,  the  double  Tafel 
slope  from  model  prediction  was  found  to  occur  around  the  same 
current  density  and  cell  voltage  as  those  measured.  The  results 
of  this  two-stage  validation  study  indicate  that  the  model  devel¬ 
oped  in  the  present  study  can  be  used  with  confidence  to  study 
the  effect  of  cathode  composition  and  structure  on  fuel  cell  perfor¬ 
mance. 

3.2.  Effect  of  cathode  composition 

Following  the  model  validation  study,  the  effect  of  cathode  com¬ 
position  on  fuel  cell  performance  was  studied.  To  this  purpose,  the 
following  four  parameters  were  varied  independently  while  hold¬ 
ing  the  others  constant: 

(1)  polymer  (Nafion)  coating  thickness  around  agglomerate 
nucleus,  8 ; 

(2)  volumetric  fraction  of  polymer  (Nafion)  in  the  agglomerate 
nucleus,  £agg; 


(3)  platinum  mass  loading,  mPt; 

(4)  platinum-carbon  mass  ratio  in  catalyst  layer  ink,  Pt|C. 

In  these  studies,  the  cathode  layer  thickness  was  held  con¬ 
stant  at  L=15  |jim  and  the  agglomerate  radius  was  held  constant 
at  ragg  =  1000  nm.  Other  baseline  values  for  the  various  parameters 
used  in  the  study  are  presented  in  Table  2.  As  the  effect  of  each 
parameter  is  analyzed,  the  optimum  value  for  that  parameter  is 
used  in  subsequent  studies,  unless  otherwise  stated. 

When  the  ionomer  (Nafion)  loading  is  changed,  it  can  result  in 
simultaneous  change  in  the  volume  fraction  of  Nafion  within  the 
agglomerate  (£agg)  and  the  ionomer  coating  thickness  (5).  FIow- 
ever,  since  there  are  infinite  many  combinations  of  the  two  model 
parameters  than  can  result  in  a  certain  prescribed  Nafion  loading 
value,  the  problem  has  been  posed  in  the  reverse  manner,  in  which 
£agg  and  8  are  varied  independently. 

3.2.1  Effect  of  polymer  coating  thickness 

Fig.  7(a)  shows  the  polarization  curves  for  L  =  15[xm, 
ragg  =  1000nm,  £agg  =  0.4,  rnPt  =  0.4mgcnrr1 2,  and  Pt|C  =  0.28  with 
the  polymer  coating  thickness  (5)  varying  from  0  to  100  nm.  The 
variation  of  the  three  volume  fractions,  namely  the  volume  fraction 
of  the  solid  (carbon  and  platinum),  the  volume  fraction  of  the 
ionomer  (Nafion),  and  the  volume  fraction  of  the  gas  pores,  is 


Fig.  7.  Effect  of  polymer  coating  thickness  ( 8 )  on  fuel  cell  performance  for  1  =  15  |jim, 
ragg  =  1000  nm,£agg  =  0.4,  mPt  =  0.4mgcm_2,  Pt|C  =  0.28.  (a)  Polarization  curves  and 
(b)  volume  fractions  of  platinum-carbon,  Nafion  and  porosity. 
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shown  separately  in  Fig.  7(b).  It  is  clear  from  the  figure  that  the  per¬ 
formance  of  the  fuel  cell  is  poor  when  the  coating  thickness  is  zero. 
As  discussed  earlier,  the  absence  of  a  coating  indicates  minimal  or 
no  contact  between  agglomerates,  and  results  in  few  continuous 
paths  within  the  catalyst  layer  for  the  protons  to  migrate.  Thus,  the 
protonic  conductivity  of  the  cathode  catalyst  layer  is  very  small 
(Eq.  (37)  and  Fig.  5)  if  5  =  0.  At  the  other  end  of  the  spectrum,  if  the 
coating  thickness  is  too  large,  the  volume  fraction  of  Nafion  within 
the  catalyst  layer  (£agg)  is  substantial  (Fig.  7(b)),  and  the  porosity 
£Cat  is  low.  This  results  in  increased  mass  transport  losses  at  the 
catalyst  layer  scale  (i.e.,  global  mass  transport  losses),  as  indicated 
by  the  downward  bend  of  the  polarization  curve  for  5  =  100  nm  at 
high  current  densities.  An  optimum  coating  thickness  is  observed 
at  higher  current  densities,  while  the  optimum  thickness  is  not  as 
discernible  at  intermediate  current  densities,  indicating  that  local 
mass  transport  limitations  due  to  increased  coating  thickness  is 
negligible.  This  observation  is  also  consistent  with  Fig.  3(d),  which 
showed  earlier  that  the  film  effectiveness  factor,  Eff2,  is  almost 
unity  irrespective  of  the  film  coating  thickness.  Based  on  these 
findings,  an  optimum  value  of  8  =  60  nm  is  used  in  all  subsequent 
studies. 

3.2.2.  Effect  of  polymer  volume  fraction  in  agglomerate  nucleus 

Fig.  8(a)  shows  the  polarization  curves  for  L=15pJm, 
ragg  =  1000nm,  8  =  60  nm,  mPt  =  0.4  mg  cm-2,  and  Pt|C  =  0.28 
with  the  polymer  (Nafion)  volumetric  fraction  in  the  agglomerate 
nucleus  (£agg)  varying  from  0.2  to  0.5.  The  variation  of  the  three 
volume  fractions,  namely  the  volume  fraction  of  the  solid  (carbon 
and  platinum),  the  volume  fraction  of  the  ionomer  (Nafion), 
and  the  volume  fraction  of  the  gas  pores,  is  shown  separately  in 
Fig.  8(b).  The  polarization  curves  show  that  at  low  to  medium 
current  densities,  the  performance  of  the  fuel  cell  improves  with 
increasing  values  of  £agg.  This  is  due  to  the  fact  that  the  polymer 
content  in  the  cathode  increases  as  £agg  increases,  thereby  reducing 
the  mass  transport  resistance  within  the  agglomerate  nucleus 
(Eq.  (4)  and  Fig.  3(a)-(c)).  At  medium  current  densities,  there  is 
no  clear  optimal  performance  for  varying  polymer  volumetric 
fraction  within  the  agglomerate  nucleus,  indicating  that  at  the 
local  or  agglomerate-scale,  increasing  the  amount  of  Nafion  within 
the  agglomerate  is  always  beneficial.  These  results  also  imply 
that  experimentally  observed  optimum  performance  with  varying 
Nafion  (polymer)  content  [3,5]  in  a  cathode  is  due  to  variation  in 
the  polymer  coating  thickness  rather  than  due  to  the  variation  of 
the  polymer  volumetric  fraction  within  the  agglomerate  nucleus. 
At  high  current  densities,  however,  global  mass  transport  losses 
increase  due  to  decrease  in  the  porosity  (£cat)  of  the  cathode  cata¬ 
lyst  layer,  as  shown  in  Fig.  8(b).  Based  on  the  results  of  this  study, 
an  optimal  value  of  £agg  =  0.4  is  used  in  all  subsequent  studies.  It  is 
worth  mentioning  at  this  point  that  the  optimum  values  chosen 
for  8  (preceding  sub-section)  and  £agg,  i.e.,  5  =  60  nm,  and  £agg  =  0.4, 
both  yield  overall  catalyst  layer  porosity  values,  £cat,  of  approxi¬ 
mately  0.28,  as  seen  in  Figs.  7(b)  and  8(b).  Incidentally,  a  value  of 
0.28  has  been  widely  cited  in  past  studies,  both  experimental  and 
numerical  [30,31,46,47,49]. 

3.2.3.  Effect  of  platinum  mass  loading 

Fig.  9(a)  shows  the  polarization  curves  for  L=15pJm, 
ragg  =  1000nm,  5  =  60  nm,  £agg  =  0.4,  and  Pt|C  =  0.28  with  plat¬ 
inum  mass  loading  (mpt)  varying  from  0.2  to  0.5  mg  cm-2.  The 
platinum  mass  loading  is  defined  as  the  mass  of  platinum  used 
per  unit  area  of  the  bi-polar  plate  through  which  the  current  is 
collected.  The  variation  of  the  three  volume  fractions,  namely 
£cat»  £n*  and  £s  is  shown  separately  in  Fig.  9(b).  At  low  to  inter¬ 
mediate  current  densities,  increase  in  platinum  mass  loading 
results  in  better  performance,  as  is  evident  from  Fig.  9(a).  This 


Fig.  8.  Effect  of  polymer  volumetric  fraction  in  agglomerate  nucleus  (£agg)  on 
fuel  cell  performance  for  L  =  15  pan,  ragg  =  1000nm,  <5  =  60nm,  mPt  =  0.4mgcirr2, 
Pt|C  =  0.28.  (a)  Polarization  curves  and  (b)  volume  fractions  of  platinum-carbon, 
Nafion  and  porosity. 

is  due  to  the  fact  that,  at  high  platinum  loading,  the  amount  of 
catalyst  surface  available  for  the  oxygen  reduction  reaction  to  take 
place  is  high  (Eq.  (9)).  At  high  current  densities,  however,  global 
mass  transport  effects  begin  to  influence  fuel  cell  performance. 
This  is  due  to  the  sharp  decrease  in  cathode  porosity,  £cat,  with 
increasing  platinum  loading,  as  shown  in  Fig.  9(b),  which  in  turn 
results  in  increased  mass  transport  resistance  at  the  cathode 
catalyst  layer  scale.  Thus,  aside  of  cost  considerations,  there  is 
a  valid  physics-based  reason  not  to  use  excessive  amounts  of 
platinum  in  PEMFC  cathodes.  Based  on  the  results  of  this  study, 
an  optimum  value  of  mPt  =  0.45  mg  cm-2  is  used  in  all  subsequent 
studies. 

3.2.4.  Effect  of  platinum-carbon  mass  ratio 

Fig.  10(a)  shows  the  polarization  curves  for  L=  15  pun, 
ragg  =  1000nm,  5  =  60  nm,  £agg  =  0.4,  and  mPt  =  0.45  mg  cm-2  with 
platinum-carbon  mass  ratio  (Pt|C)  varying  from  0.25  to  0.4.  Chang¬ 
ing  the  Pt|C  ratio  with  a  fixed  platinum  mass  loading  implies  that 
the  amount  of  carbon  within  the  agglomerates  is  changed.  Further¬ 
more,  since  the  volume  fraction  of  Nafion  within  the  agglomerate, 
£agg,  is  also  to  be  kept  unchanged,  change  in  the  amount  of  car¬ 
bon  must  also  be  accompanied  by  a  changing  in  the  amount  of 
Nafion  within  the  catalyst  layer.  The  variation  of  the  three  relevant 
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Fig.  9.  Effect  of  platinum  mass  loading  (mPt)  on  fuel  cell  performance  for  I  =  15  p,m, 
ragg  =  1000  nm,  5  =  60  nm,  £agg  =  0.4,  Pt|C  =  0.28.  (a)  Polarization  curves  and  (b)  volume 
fractions  of  platinum-carbon,  Nation  and  porosity. 


global  volume  fractions  is  shown  in  Fig.  10(b).  It  can  be  seen  that 
at  low  platinum-carbon  mass  ratios,  the  porosity  of  the  cathode 
is  low.  Based  on  the  above  discussion,  low  Pt|C  ratio  implies  large 
amounts  of  carbon,  which  also  implies  large  amounts  of  Nation. 
Thus,  low  Pt|C  ratio  produces  poor  performance  due  to  increased 
mass  transport  resistance  to  diffusion  of  gases  through  the  pores,  as 
indicated  by  the  strong  bending  of  the  dotted  red  curve  at  high  cur¬ 
rent  densities.  At  the  other  end,  excessively  high  values  of  Pt|Cratio 
imply  that  there  is  too  little  Nation  within  the  cathode,  thereby 
reducing  the  protonic  conductivity  of  the  cathode.  Thus,  we  see 
that  higher  Pt|C  values  give  rise  to  lower  performance,  irrespec¬ 
tive  of  the  current  density.  Based  on  the  results  of  the  current 
study,  an  optimum  value  of  Pt|C  =  0.28  is  used  in  all  subsequent 
studies. 

3.3.  Effect  of  cathode  structure 

In  order  to  study  the  effect  of  cathode  structure  on  fuel  cell 
performance,  the  following  two  parameters  were  varied  indepen¬ 
dently: 

(1)  agglomerate  radius  (ragg); 

(2)  cathode  layer  thickness  (L). 


3.3 A.  Effect  of  agglomerate  radius 

Fig.  11(a)  shows  the  polarization  curves  for  L  =  15  [xm,  5  =  60  nm, 
£agg  =  0.4,  /7ipt  =  0.45  mg  cm-2,  and  Pt|C  =  0.28  with  agglomerate 
radius  (ragg)  varying  from  600  to  1200  nm.  In  this  first  study,  the 
coating  thickness  is  kept  unchanged  even  though  the  agglomerate 
radius  is  changed.  This  implies,  that  the  volume  fraction  of  Nation 
within  the  catalyst  layer,  £N,  will  change,  as  will  the  overall  poros¬ 
ity  of  the  cathode  catalyst  layer.  The  variation  of  the  three  relevant 
global  volume  fractions  is  shown  in  Fig.  11(b).  Fig.  11(a)  shows  that 
as  the  agglomerate  radius  increases,  the  performance  of  the  fuel  cell 
deteriorates.  Essentially,  smaller  agglomerates  pose  smaller  resis¬ 
tance  to  mass  transport  at  the  agglomerate  (local)  scale,  as  was 
also  illustrated  in  Fig.  3(a)-(c).  Thus,  the  performance  monotoni- 
cally  improves  with  decrease  in  agglomerate  radius,  as  predicted 
by  Sun  et  al.  [18]  and  Secanell  et  al.  [19].  Flowever,  experimental 
results  [27]  suggest  that  performance  deteriorates  if  the  agglomer¬ 
ate  size  is  too  small.  Some  authors  [15,16,27]  cite  this  as  Knudsen 
effects.  They  argue  that  small  agglomerate  sizes  also  result  in  very 
small  pores,  within  which  Knudsen  effects  are  manifested.  How¬ 
ever,  since  ballistic  transport  poses  lower  resistance  than  diffusive 
transport,  it  is  not  clear  if  Knudsen  effects  can  account  for  this 
deterioration  in  performance,  i.e.,  increase  in  mass  transport  resis¬ 
tance.  It  is  our  contention  that  this  deterioration  of  performance  for 


Fig.  10.  Effect  of  platinum-carbon  mass  ratio  in  catalyst  layer  ink  (Pt|C)  on  fuel  cell 
performance  forL=15pim,  ragg  =  1000nm,  5  =  60 nm,  eagg  =  0.4,  mPt  =  0.45 mgcm-2. 
(a)  Polarization  curves  and  (b)  volume  fractions  of  platinum-carbon,  Nation  and 
porosity. 
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Fig.  11.  Effect  of  agglomerate  radius  (ragg)  on  fuel  cell  performance  for  L  =  15  fxm, 
<5  =  60nm,  £agg  =  0.4,  mpt  =  0.45mgcm-2,  Pt|C  =  0.28.  (a)  Polarization  curves  and  (b) 
volume  fractions  of  platinum-carbon,  Nation  and  porosity. 


small  agglomerate  sizes  is  due  to  increased  tortuosity  of  the  pores. 
Experimental  and  theoretical  studies  [28]  clearly  indicate  that  a  box 
packed  randomly  with  large  rigid  spheres  pose  smaller  resistance 
to  diffusion  than  a  box  packed  with  small  spheres  with  the  same 
overall  porosity.  The  model  used  here  (Eq.  (23)  and  Fig.  4)  clearly 
captures  this  physics.  As  a  result,  our  model  predictions  indeed 
show  (Fig.  11(a))  that  at  high  current  densities,  smaller  sizes  do 
not  result  in  improved  performance.  The  optimum  is  governed  by 
the  balance  between  decreased  local  mass  transport  resistance  due 
to  smaller  agglomerate  size  and  an  increased  global  mass  transport 
resistance  due  to  increased  tortuosity. 

In  the  first  set  of  studies,  just  described,  change  in  the  agglom¬ 
erate  radius  also  resulted  in  change  of  the  three  relevant  global 
volume  fractions.  In  other  words,  the  structural  parameter  affected 
the  compositional  parameters.  In  order  to  isolate  the  effect  of 
the  structural  parameter  on  fuel  cell  performance,  agglomerate 
radius  needs  to  be  varied  in  such  a  way  that  the  volume  frac¬ 
tions  of  platinum-carbon,  Nation  and  the  porosity  of  the  cathode 
remain  constant.  This  is  achieved  by  varying  the  agglomerate 
radius  while  keeping  the  ratio  between  polymer  coating  thick¬ 
ness  and  the  agglomerate  radius  constant  (see  Fig.  12(b)).  Fig.  12(a) 
shows  the  polarization  curves  for  L  =  15  [Jim,  8/ragg  =  0.06,  sagg  =  0.4, 
mPt  =  0.45  mg  cm-2,  and  Pt|C  =  0.28  with  agglomerate  radius  (ragg) 
varying  from  50  nm  to  1000  nm.  From  Fig.  12(a),  it  is  clear  that  while 


performance  improves  with  decreasing  agglomerate  radius,  at  very 
low  values  of  ragg,  there  is  no  significant  change  in  performance,  as 
explained  in  the  preceding  paragraph.  In  fact,  at  high  current  den¬ 
sities,  the  performance  for  ragg  =  100nm  is  seen  to  be  marginally 
better  than  that  for  ragg  =  50  nm.  This  implies  that  cathode  tortuos¬ 
ity  has  a  significant  influence  on  fuel  cell  performance  at  low  values 
of  agglomerate  radius. 

3.3.2.  Effect  of  cathode  catalyst  layer  thickness 

Fig.  13(a)  shows  the  polarization  curves  for  ragg  =  200nm, 
8/ragg  =  0.06,  £agg  =  0.4,  and  Pt|C  =  0.28  for  cathode  layer  thicknesses 
of  L  =  7.5,  15,  and  30  [xm.  The  cathode  catalyst  layer  thickness  is  a 
global  structural  parameter,  and  it  does  not  affect  the  three  rele¬ 
vant  porosities,  as  is  clear  from  Fig.  13(b).  It  can  be  seen  that  the 
performance  is  superior  for  L  =  15  |xm  compared  to  the  L  =  7.5  |xm 
case.  This  is  due  to  the  fact  that  the  reactant  gases  do  not  have  as 
many  reaction  sites  in  the  latter  case  as  in  the  former  case.  Com¬ 
paring  L  =  15  and  30  pan,  it  can  be  seen  that  the  latter  thickness 
provides  better  performance  at  low  and  intermediate  current  den¬ 
sities.  Due  to  increased  mass  transport  resistance  for  L  =  30  pan,  at 
higher  current  densities,  the  reactant  gases  are  unable  to  reach 
all  the  reaction  sites,  and  the  catalyst  at  the  inner  half  of  the 
catalyst  layer  (i.e.,  adjacent  to  the  membrane)  largely  remain  under¬ 
utilized. 


Current  Density  A  (cm  2) 


Fig.  12.  Effect  of  agglomerate  radius  (ragg)  on  fuel  cell  performance  for  1  =  15  |jim, 
<5/ragg  =  0.06,  £agg  =  0.4,  mPt  =  0.45  mg  cm-2,  Pt|C  =  0.28.  (a)  Polarization  curves  and  (b) 
volume  fractions  of  platinum-carbon,  Nafion  and  porosity. 
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Fig.  13.  Effect  of  cathode  catalyst  layer  thickness  (I)  on  fuel  cell  performance  for 
rag g  =  200  nm,  <5/ragg  =  0.06,  eagg  =  0.4,  Pt|C  =  0.28.  (a)  Polarization  curves  and  (b)  vol¬ 
ume  fractions  of  platinum-carbon,  Nation  and  porosity. 

4.  Summary  and  conclusions 

With  the  advent  of  next-generation  techniques  for  micro-/nano- 
fabrication,  it  will  be  possible  to  engineer  fuel  cell  catalyst  layers 
(electrodes)  to  desired  compositions  and  structures.  If  a  model 
is  available  that  is  general  enough  to  answer  the  question,  even 
qualitatively,  as  to  what  structure  and  composition  is  best  for  opti¬ 
mum  performance,  its  impact  on  the  advancement  of  fuel  cell 
technology  could  be  unprecedented.  The  model  developed  in  this 
study  fulfills  this  critical  need.  The  sub-grid  scale  cathode  cata¬ 
lyst  model,  based  on  the  flooded  agglomerate  concept,  was  first 
validated  against  previously  published  experimental  data.  Follow¬ 
ing  the  validation  phase,  the  effects  of  Nation  loading,  platinum 
loading,  platinum-carbon  ratio,  agglomerate  size  and  cathode  layer 
thickness  on  the  performance  of  a  polymer  electrolyte  membrane 
fuel  cell  were  studied.  Based  on  these  studies,  the  following  con¬ 
clusions  may  be  drawn: 

(1)  Variation  of  engineering  parameters  such  as  platinum  loading 
or  Nation  loading  leads  to  simultaneous  change  in  several  of  the 
model  parameters— most  notably  the  volume  fractions  of  the 
pores,  Nation  and  solid.  Thus,  the  performance  of  the  overall 
fuel  cell  is  a  manifestation  of  a  complex  interaction  of  several 
model  parameters. 


(2 )  The  effect  of  each  parameter  on  the  overall  performance  of  a  fuel 
cell  is  different  at  different  operating  points  (current  densities). 
While  local  (or  agglomerate-scale)  mass  transport  effects  can 
be  understood  simply  by  carefully  studying  the  effectiveness 
factors  shown  in  Fig.  3,  global  mass  transport  losses  require 
incorporation  of  the  sub-grid  scale  catalyst  layer  model  into 
full-scale  CFD  calculations. 

(3)  At  medium  current  densities,  variation  of  compositional 
parameters  usually  affects  performance  due  to  local  mass  trans¬ 
port  and  electrical  conduction  limitations,  while  at  high  current 
densities,  performance  is  primarily  affected  due  to  global  mass 
transport  limitations  introduced  arising  from  varying  porosity. 

(4)  The  presence  of  an  optimal  performance  with  varying  Nation 
content  in  the  cathode  is  more  due  to  the  local  agglomerate- 
level  mass  transport  and  conductivity  losses  in  the  polymer 
coating  around  the  agglomerates  than  due  to  the  amount  of 
Nation  within  the  agglomerate. 

(5)  In  addition  to  cost  constraints,  platinum  mass  loading  needs 
to  be  at  a  moderate  level  in  order  to  optimize  fuel  cell  perfor¬ 
mance. 

(6)  Cathode  tortuosity  has  a  significant  effect  on  fuel  cell  per¬ 
formance  at  low  values  of  agglomerate  radius.  Thus,  it  is  not 
desirable  to  design  cathodes  in  which  the  carbon  particle  sizes 
(or  pore  sizes)  are  too  small.  While  this  may  result  in  reduced 
agglomerate-scale  mass  transport  losses,  it  is  undesirable  due 
to  poor  diffusion  of  reactants  through  the  pores. 

In  the  future,  the  study  will  be  extended  to  include  the  effect  of 

liquid  water.  Flowever,  the  authors  hope  that  this  study  has  been 

able  to  lay  a  solid  foundation  for  further  advancement  in  this  area. 
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